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Abstract 



A general theorem is rigorously proved for the case, when an observable is 
a sum of linearly independent terms: The dispersion of a global observable is 
normal if and only if all partial dispersions of its terms are normal, and it is 
anomalous if and only if at least one of the partial dispersions is anomalous. 
This theorem, in particular, rules out the possibility that in a stable system 
with Bose-Einstein condensate some fluctuations of either condensed or non- 
condensed particles could be anomalous. The conclusion is valid for arbitrary 
systems, whether uniform or nonuniform, interacting weakly or strongly. The 
origin of fictitious fluctuation anomalies, arising in some calculations, is elu- 
cidated. 
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The problem of fluctuations of observable quantities is among the most important ques- 
tions in statistical mechanics, being related to the fundamentals of the latter. A great 
revival, in recent years, of interest to this problem is caused by intensive experimental and 
theoretical studies of Bose-Einstein condensation in dilute atomic gases (see e.g. reviews [1- 
3] ) . The fact that the ideal uniform Bose gas possesses anomalously large number-of-particle 
fluctuations has been known long ago [4,5], which has not been of much surprise, since such 
an ideal gas is an unrealistic and unstable system. However, as has been recently suggested 
in many papers, similar anomalous fluctuations could appear in real interacting Bose sys- 
tems. A number of recent publications has addressed the problem of fluctuations in Bose 
gas, proclaiming controversial statements of either the existence or absence of anomalous 
fluctuations (see discussion in review [6]). So that the issue has not been flnally resolved. 
In the present paper, the problem of fluctuations is considered from the general point of 
view, independent of particular models or calculational methods. A general theorem is rig- 
orously proved, from which it follows that there are no anomalous fluctuations in any stable 
equilibrium systems. 

It is worth stressing that no phase transitions are considered in this paper. As can be 
easily inferred from any textbook on thermodynamics or statistical mechanics, the points 
of phase transitions are, by deflnition, the points of instability. A phase transition occurs 
exactly because one phase becomes unstable and has to change to another stable phase. It 
is well known that at the points of second-order phase transitions fluctuations do become 
anomalous, yielding divergent susceptibilities, as it should be at the points of instability. 
After a phase transition has occurred, the system, as is also well known, becomes stable 
and susceptibilities go finite. However, in many papers on Bose systems, the claims are 
made that fiuctuations remain anomalous far below the condensation point, in the whole 
region of the Bose-condensed system. As is shown below, these claims are incorrect, since 
such a system with anomalous fluctuations possesses a divergent compressibility, thus, being 
unstable. 

Observable quantities are represented by Hermitian operators from the algebra of ob- 
servables. Let A be an operator from this algebra. Fluctuations of the related observable 
quantity are quantifled by the dispersion 

A\A) = < > - < i >2 , (1) 

where < . . . > implies equilibrium statistical averaging. The dispersion itself can be treated 
as an observable quantity, which is the average of an operator {A— < A >)^, since each 
dispersion is directly linked to a measurable quantity. For instance, the dispersion for the 
number-of-particle operator deflnes the isothermal compressibility 

in which P is pressure, p = N/V is density, N =< N >, V is volume, and T temperature. 
The dispersion A^(iV) is also connected with the sound velocity s through the equation 
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where m is particle mass, P pressure, and with the central value 



5(0)=M.T..= ^4 = ^ (4) 

of the structural factor 

S{k) -^I + Pj [9{r) - l]e-^'^-^c/r , 

where g{r) is the pair correlation function. The fluctuations of the Hamiltonian H charac- 
terize the specific heat 




N\dT)y NkeT^ ' 

where E =< H > is internal energy. In magnetic systems, with the Zeeman interaction 
~A*o Z^i B • Sj, the longitudinal susceptibility 

= 1 = ^'(^«) 

is related to the fluctuations of the magnetization Ma =< Ma >, where Mq = //q J2iLi S^. 

All relations (2) to (6) are exact and hold true for any equilibrium system. The stability 
conditions for such systems require that at any finite temperature, except the points of phase 
transitions, quantities (2) to (6) be positive and finite for all N, including the thermodynamic 
limit, when N ^ oo. At the phase transition points, these quantities can of course be 
divergent, since, as is well known, the phase transition points are the points of instability. 
Summarizing this, we may write the general form of the necessary stability condition as 

0<^<oc. (7) 

which must hold for any stable equilibrium systems at finite temperature and for all N, 
including the limit N ^ oo. The value A'^{A)/N can become zero only at T = 0. Condition 
(7) is nothing but a representation of the well known fact that the susceptibilities in stable 
systems are positive and finite. 

The stability condition (7) shows that the dispersion A^(A) has to be of order N. When 
A^(A) ~ N, one says that the dispersion is normal and the fiuctuations of an observable A 
are normal, since then the stabihty condition (7) is preserved. But when A'^{A) ~ A'"'*, with 
a > 1, then such a dispersion is called anomalous and the fluctuations of A are anomalous, 
since then A'^{A)/N ~ A^"~^ ^ cxo as A^ oo, hence the stability condition (7) becomes 
broken. A system with anomalous fiuctuations is unstable. For example, an ideal uniform 
Bose gas, with A^(iV) ~ A^^, is unstable [4-6]. 

Thus, in any stable equilibrium system, the fluctuations of global observables must be 
normal, A^(A) ~ A^. The situation becomes more involved, if the operator of an observable 
is represented by a sum 

i = ^i, (8) 
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of Hcrmitian terms Ai, and one is interested in the fluctuations of tlie latter. Then the 
intriguing question is: Could some partial dispersions A^(ylj) be anomalous, while the total 
dispersion A^(A) remaining normal? Exactly such a situation concerns the systems with 
Bose condensate. Then the total number of particles = A^o + is a sum of the numbers of 
condensed, A^o, and noncondensed, A^i, particles. We know that for a stable system A^(7V) 
must be normal. But could it happen that at the same time either A^(7Vo) or A^(7Vi), or 
both would be anomalous, as is claimed by many authors? 

Considering the sum (8), it is meaningful to keep in mind a nontrivial case, when all Ai 
are linearly independent. In the opposite case of linearly dependent terms, one could simply 
express one of them through the others and reduce the number of terms in sum (8). The 
consideration also trivilizes if some of are c-numbers, since then A^(c) =0. 

The dispersion of operator (8) reads as 

A^A) = J2A'{Ai) + 2j2coviA,A,) , (9) 

i i<j 

where the covariance 

is introduced. The latter is symmetric, cov{Ai, Aj) = cov{Aj, Ai). The dispersions are, 
by definition, non-negative, but the covariances can be positive as well as negative. One 
might think that an anomalous partial dispersion A^(Ai) could be compensated by some 
covariances, so that the total dispersion A'^{A) would remain normal. This is just the way of 
thinking when one finds an anomalous dispersion of condensed, A^(iVo), or noncondensed, 
A^(A^i), particles, presuming that the system as a whole could remain stable, with the 
normal total dispersion A^(A^). However, the following theorem rules out such hopes. 

Theorem. The total dispersion (9) of an operator (8), composed of linearly independent 
Hermitian operators, is anomalous if and only if at least one of the partial dispersions is 
anomalous, with the power of the total dispersion defined by that of its largest partial 
dispersion. Conversely, the total dispersion is normal if and only if all partial dispersions 
A^(Ai) are normal. 

Proof. First of all, we notice that it is sufficient to prove the theorem for the sum of 
two operators, for which 

A\Ai + Aj) = A\Ai) + A\Aj) + 2 cov(A, Aj) , (10) 

where i ^ j. This is because any sum of terms more than two can always be represented as a 
sum of two new terms. Also, we assume that both operators in Eq. (10) are really operators 
but not c-numbers, since if at least one of them, say Aj = c, is a c-number, then Eq. (10) 
reduces to a simple equality A'^{Ai + c) = A'^{Ai) of positive (or semipositive) quantities, 
both of which simultaneously are either normal or anomalous. 
Introduce the notation 

aij = cov(Aj, Aj) . (11) 



= - < AiAj + AjAi > - <Ai>< Aj > 
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As is evident, an = A^(Aj) > and aij = aji. The set of elements aij forms the covariance 
matrix [cij], which is a symmetric matrix. For a set of arbitrary real- valued numbers, Xi, 
with i = 1,2, ... ,n, where n is an integer, one has 



n "1 2 

< V (Ai- <Ai>) Xi 



> = ^aijXiXj > . (12) 
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The right-hand side of equality (12) is a semipositive quadratic form. From the theory of 
quadratic forms [7] one knows that a quadratic form is semipositive if and only if all principal 
minors of its coefficient matrix are non-negative. Thus, the sequential principal minors 
of the covariance matrix [cTij], with i,j — 1,2, ... ,n, are all non- negative. In particular, 
o'iiO'jj — <^ij<^ji > 0. This, owing to the symmetry aij — aji, transforms to the inequality 
f^tj < CTnajj. Then the correlation coefficient 

Xij = aij/y/anajj (13) 

possesses the property Xfj < 1. 

The equality X^j = 1 holds true if and only if Ai and Aj are linearly dependent. The 

sufficient condition is straightforward, since if Aj = a + bAi, where a and b are any real 
numbers, then aij — ban and ajj — b'^an, hence Xij — b/\b\, from where A?- = 1. To prove 
the necessary condition, assume that A? = 1. This implies that Xij = ±1. Consider the 
dispersion 

'"^ ^' ±^U2(1±A,,)>0. 



\ y^ii \^33 / 

The value Aj^ = 1 is possible then and only then, when 

. o I Ai 




The dispersion can be zero if and only if 



Ai Aj 



'^33 



const , 



that is, the operators Aj and Aj are linearly dependent. Similarly, the value Xij — —1 is 
possible if and only if 

— — -I — const , 



'^33 



which again means the linear dependence of the operators Ai and Aj. As far as these 
operators are assumed to be linearly independent, one has A^^ < 1. The latter inequality is 
equivalent to aij < anajj, which, in agreement with notation (11), gives 

\cov{A, Aj)\' <A\Ai)A\Aj) . (14) 
The main relation (10) can be written as 
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(15) 



where, as is shown above, \Xij\ < 1. Altogether there can occur no more than four following 
cases. First, when both partial dispersions an = A^(yij) and ajj = A^(Aj) are normal, so 
that an ~ N and ajj ~ N. Then from Eq. (15) it is evident that the total dispersion 
A.^{Ai + Aj) ~ is also normal. Second, one of the partial dispersions, say an ~ A^, is 
normal, but another is anomalous, ajj ~ A^", with a > 1. Prom Eq. (15), because of 
(1 + a)/2 < a, one has A^(Aj + Aj) ~ A^", so that the total dispersion is anomalous, having 
the same power a as ajj. Third, both partial dispersions are anomalous, an ~ A^"^ and 
ajj ~ N""^, with different powers, say 1 < ai < a2. Prom Eq. (15), taking into account that 
(cKi + q;2)/2 < q;2, we get A.'^(Ai + Aj) ~ A^"^^ ^^^^ jg^ ^j^g total dispersion is also anomalous, 
with the same power a2 as the largest partial dispersion ajj. Fourth, both partial dispersions 
are anomalous, an = cjN" and ajj — CjN"', where Cj > and Cj > 0, with the same power 
a. Then Eq. (15) yields 

A^{Ai + Aj) = CijN'' 

with 

Cij = {ci - Cj)"^ + 2ciCj{l + Xij) > , 

which is strictly positive in view of the inequality |Ajj| < 1. Hence, the total dispersion is 
anomalous, with the same power a as both partial dispersions. After listing all admissible 
cases, we see that the total dispersion A^(Aj + Aj) is anomalous if and only if at least one 
of the partial dispersions is anomalous, with the power of A^ of the total dispersion being 
equal to the largest power of partial dispersions. Oppositely, the total dispersion is normal 
if and only if all partial dispersions are normal. This concludes the proof of the theorem. 

As an example, let us consider a Bose system with Bose-Einstein condensate, whose 
total numbcr-of-particlc operator A^ = Aq + iVi consists of two terms, corresponding to 
condensed, A'o; and noncondensed, Ni, particles. Since for a stable system the dispersion 
A^(7V) is normal, then from the above theorem it follows that both dispersions A^(7Vo) as 
well as A^(A^i) must be normal. No anomalous fluctuations can exist in a stable system, 
neither for condensed nor for noncondensed particles. This concerns any type of stable 
systems, either uniform or nonuniform. And this result does not depend on the method of 
calculations, provided the latter are correct. 

How then could one explain the appearance of numerous papers claiming the existence of 
anomalous fluctuations in Bose-condensed systems in the whole region far below the critical 
point? If these anomalous fluctuations would really exist, then the compressibility would be 
divergent everywhere below the critical point. A system, whose compressibility is divergent 
everywhere in its region of existence, as is known from any textbook on statistical mechanics, 
is unstable. Such anomalous fluctuations are usually obtained as follows. One considers a 
low-temperature dilute Bose gas, at T <^ Tc, when the Bogolubov theory [8] is applicable. 
In the frame of this theory, one calculates the dispersion A^(A^i) of noncondensed particles, 
where Ni = J2kjto (^l^k- To find < A"^ >, one needs to work out the four-operator expression 
< alaka^^aq >, or after employing the Bogolubov canonical transformation = Ukbk + 
"^-fc^-fc! ^'^ consider < bl.bkb]^bg >. Such four-operator expressions arc treated by invoking 
the Wick decoupling. Then one finds that the dispersion A^(Ai) diverges a.s N J dk/k'^. 
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Discretizing the phonon spectrum, one gets A^(7Vi) ~ N^/^. Another way [6] could be 
by hmiting the integration by the minimal kmin = '^/L, where L ~ N^'''^ is the system 
length. Then again A^(iV'i) ~ N^^^. In any case, one obtains anomalous fluctuations 
of noncondensed particles. This result holds true for both canonical and grand canonical 
ensembles, which is a direct consequence of the Bogolubov theory [8] . But the anomalous 
behaviour of A^(iVi), according to the theorem proved above, immediately leads to the same 
anomalous behaviour of the total dispersion A^(A^), which would imply the instability of the 
system, since then compressibility (2) and structural factor (4) become divergent. As far as 
the system is assumed to be stable, there should be something wrong in such calculations. 

The drawback of these calculations is in the following. One of the basic points of the 
Bogolubov theory is in omitting in the Hamiltonian all terms of orders higher than two with 
respect to the operators at of noncondensed particles. This is a second- order theory with 
respect to at- It means that the same procedure of keeping only the terms of second order, 
but ignoring all higher-order terms, must be done in calculating any physical quantities. In 
considering < TVf >, one meets the fourth-order terms with respect to Ofe. Such fourth-order 
terms are not defined in the Bogolubov theory. The calculation of the fourth-order products 
in the second-order theory is not self-consistent. This inconsistency leads to incorrect results. 

A correct calculations of A^(A^) in the frame of the Bogolubov theory can be done in the 
following way [9]. Writing down the pair correlation function, one should retain there only 
the terms not higher than of the second order with respect to a^, omitting all higher-order 
terms. Then for a uniform system one gets 

^(r) = 1 + - /" (< aluk > + < aka_k >) e*"" . 



This gives us the structural factor S{0) = ksT/mc^, c = y (p/m)$o, $o = / $(r)(ir, where 

$(r) is an interaction potential. Because of the exact relation A^(iV) = NS{0), we find the 
dispersion A^(iV) = NksT /mc^, which is, of course, normal, as it should be for a stable 
system. And since the total dispersion A^^TV) is normal, the theorem tells us that both the 
partial dispersions, A^(A''o) as well as A^(iVi), must also be normal. Note that if in defining 
the pair correlation function g{T) one would retain the higher-order terms, one would again 
get the anomalous total dispersion A^(A^). 

It is easy to show that the same type of fictitious anomalous fluctuations appear for ar- 
bitrary systems, if one uses the second-order approximation for the Hamiltonian but intends 
to calculate fourth-order expressions. This is immediately evident from the analysis of sus- 
ceptibilities for arbitrary systems with continuous symmetry, as is done by Patashinsky and 
Pokrovsky [10] (Chapter IV), when the system Hamiltonian is restricted to hydrodynamic 
approximation. Following Ref. [10], one may consider an operator A — A{ip) being a func- 
tional of a fleld if. Let this operator be represented as a sum A — Aq + Ai^vo. which the flrst 
term is quadratic in the fleld 99, so that Aq ~ (p'^^p, while the second term depends on the 
fleld fluctuations 5ip as Ai ~ 5ip^5ip. Keeping in the Hamiltonian only the second-order fleld 
fluctuations is equivalent to the hydrodynamic approximation. The dispersion A^(yl) ~ Nx 
is proportional to a /on^i^uf/ma/ susceptibility %. The latter is given by the integral / C(r)(ir 
over the correlation function C(r) = g{v) — 1, where g{r) is a pair correlation function. 
Calculating A^(A), one meets the fourth-order term < Sip~^S(pS(p^Sip >. If this is treated 
by invoking the Wick decoupling and the quadratic hydrodynamic Hamiltonian, one gets 
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C(r) ~ l/r^(''~^) for any dimensionality d > 2. Consequently, x ~ I C{r)dr Ni^-'^)/^ for 
2 < d < A, and the dispersion A^(A) ~ Nx ~ A^("'+i)/3. For c? = 3, this results in the 
anomalous dispersion A^(yl) ~ A*"^/^. If this would be correct, this would mean, according 
to the necessary stability condition (7), that the system is unstable. That is, there could 
not exist any stable systems with continuous symmetry, such as magnetic systems or liquid 
helium. Of course, we know that such systems perfectly exist, but the above contradiction 
has arisen solely due to an inconsistent calculational procedure, when the fourth-order term 
< d(p~^Sip6(p~^6(p > was treated in the frame of the hydrodynamic approximation, which is a 
second-order theory with respect to Sip. 

Concluding, there are no anomalous fluctuations of any physical quantities in arbitrary 
stable equilibrium systems. Fictitious anomalous fluctuations in realistic systems far outside 
any points of phase transitions might appear only due to drawbacks in a calculational pro- 
cedure. The absence of anomalous fluctuations follows from the general theorem, rigorously 
proved in this paper. 

It is important to stress that not only the anomalous fluctuations as such signify the 
occurrence of instability, though they are the explicit signals of the latter. But also one 
should not forget that the dispersions for the operators of observables are directly related to 
the corresponding susceptibilities, as in Eqs. (2) to (6). It is the wrong behaviour of these 
susceptibilities, which manifests the instability. 

Thus, the anomalous number-of-particle fluctuations are characterized by the anomalous 
dispersion A^(iV). The latter is connected, through the general and exact relation (2), with 
the isothermal compressibility kt- The anomalous dispersion A^(A^) implies the divergence 
of this compressibihty. Prom the definition of the compressibihty kt = —{1/V) {dV/ dP) , it is 
obvious that its divergence means the following: An infinitesimally small positive fluctuation 
of pressure abruptly squeezes the system volume to a point. Respectively, an inflnitesimally 
small negative fluctuation of pressure suddenly expands the system volume to inflnity. It 
is more than evident that a system which is unstable with respect to infinitesimally small 
fiuctuations of pressure, immediately collapsing or blowing up, has to be termed unstable. 
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